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RESUMEN 

Presentamos espectrofotometria de alta resolucion de NGC 5461 y NGC 5471, 
dos regiones HII gigantes en la galaxia MlOl. Los datos se obtuvieron con el teles- 
copic de 2.1 m del Observatorio Astronomico Nacional en San Pedro Martir, Baja 
California. Medimos las intensidades de un conjunto de Ifneas de recombinacion de 
hidrogeno y de helio, asf como Ifneas prohibidas de un gran niimero de iones. Cal- 
culamos las condiciones fisicas en las dos nebulosas y calculamos las abundancias 
quimicas totales tomando en cuenta las abundancias ionicas observadas asf como 
las no observadas; estas ultimas las estimamos a partir de factores de correccion 
de ionizacion («c/'s). Para NGC 5461 los ic/'s se basan en un niodelo detallado 
de fotoionizacion disenado especfficamente para este objeto (Luridiana & Peimbert 
2001), mientras que para NGC 5471 se obtuvieron a partir de un modelo de fo- 
toionizacion para NGC 2363 (Luridiana, Peimbert, & Leitherer 1999), region que 
muestra un grado de ionizacion muy semejante al de NGC 5471. Los icjs asf de- 
terminados los comparamos con aquellos que se obtienen a partir de las formulas 
de Mathis & Rosa (1991). Dicha comparacion muestra importantes discrepancias 
para algunos de los elementos incluyendo nitrogeno, neon y cloro. 

ABSTRACT 

We present high resolution spectroscopic data of the two giant extragalactic 
H II regions NGC 5461 and NGC 5471 in MlOl, which have been obtained with the 
2.1-m telescope of the Observatorio Astronomico Nacional in San Pedro Martir, 
Baja California. We measured the intensities of several H and He recombination 
lines, and of forbidden lines of a large number of ions. We calculate the physical 
conditions in the two nebulae with a large number of diagnostics and determine 
their chemical abundances by applying ionization correction factors (ic/'s) to the 
observed ionic abundances. For NGC 5461, the ic/'s are based on a tailored pho- 
toionization model of the region (Luridiana & Peimbert 2001), while for NGC 5471 
they are computed from those predicted by a photoionization model of NGC 2363 
(Luridiana, Peimbert, & Leitherer 1999), a region which is similar to NGC 5471 in 
the ionization structure. For both regions, the ic/'s are compared to those com- 
puted following the prescriptions by Mathis & Rosa (1991). Such comparison shows 
large discrepancies for several elements, including nitrogen, neon, and chlorine. 

Key Words: H II REGIONS ISM: ABUNDANCES ISM: INDIVID- 
UAL (NGC 5461, NGC 5471) 



1. INTRODUCTION 

NGC 5461 and NGC 5471 are two giant extra- 
galactic HII regions (GEHRs), located in the spiral 
galaxy M 101 (NGC 5457). Due to their prominence. 
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the two GEHRs have been extensively studied by 
many authors. Low- and medium-resolution spec- 
troscopic data in several wavelength intervals, rang- 
ing from the IR to the UV, of one or both regions 
have been published by Rayo, Peimbert, & Torres- 
Peimbert (1982), McCaU, Rybski, & Shields (1985), 
Torres-Peimbert, Peimbert, & Fierro (1989), Skill- 
man & Israel (1988), Castaheda, Vflchez, & Copetti 
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TABLE 1 
JOURNAL OF OBSERVATIONS 



Date 


Spectral Range 


Orders 


Designation 


Slit Size 


Exposure Times (seconds) 
NGC 5461 NGC 5471 


96 Jun 10 
96 Jun 13 


6280 9100 A 
3450-6650 A 


25-35 
34-64 


Near IR 
Blue 


39.9 X 4" 
13.3 X 4" 


60, 3 X 1200 60, 2 x 600, 1200 
180, 2 X 900 180, 2 x 600, 900 



(1992), Rosa & Benvcnuti (1994), and Garnctt et al. 
(1999). Theoretical models covering different prop- 
erties of the two regions have been described by 
Evans (1986), Castaficda ct al. (1992), and Lurid- 
iana & Peimbert (2001). The chemical composition 
of the two regions has been determined by Torres- 
Pcimbert et al. (1989), who aimed at determining 
the chemical gradients of the parent galaxy. 

In the present paper, we present high-resolution 
spectroscopic data of both H II regions covering a 
wide spectral range, from near-ultraviolet to near- 
infrared. The use of the echelle cross disperser pro- 
duces deep high-resolution spectra of a large number 
of lines in the whole optical domain, overcoming the 
problem of line-blending, and allowing the applica- 
tion of a large number of physical-condition diag- 
nostics. The structure of the paper is as follows: in 
§ 2 we describe the observations and the reduction 
procedure; in § 3 we present the line intensities of 
the two regions; in § 4 we derive the physical condi- 
tions in the two nebulae, and calculate their chem- 
ical abundances by means of standard techniques; 
and finally, in § 5 we discuss our results and their 
implications. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were carried out with the 2.1-m 
telescope of the Obscrvatorio Astronomico Nacional 
in San Pedro Martir, Baja California, Mexico, in 
June 1996. The telescope was in its f/7.5 configu- 
ration. High resolution CCD spectra were obtained 
using the REOSC Echelle Spectrograph; see Levine 
& Chakrabarty (1994) for a description of the gen- 
eral characteristics of this instrument. The echelle 
gives a resolution of 0.234 A pixel" ^ at Ha using the 
University College of London (UCL) camera and a 
CCD-Tek chip of 1024 x 1024 pixels with a 24 /zm^ 
pixel size. The spectral resolution is 0.5 A FWHM 
and the accuracy in the wavelength determination of 
emission lines is 0.1 A. 

For both regions, we obtained spectra in two par- 
tially overlapped wavelength intervals covering the 
range from 3450 to 9100 A. Typically, three or four 



individual exposures were added to obtain the final 
spectra in each interval. Slits covering 13.3" x 4" in 
the blue and 39.9" x 4" in the near-infrared (NIR) 
were used to avoid overlapping between orders in the 
spatial direction. The slit orientation was east-west 
in all cases. A journal of the observations is pre- 
sented in Table 1. 

The atmospheric refraction depends on the ob- 
served wavelength; therefore the region of the sky 
included in the slit varies with wavelength. This 
problem is called the atmospheric differential refrac- 
tion problem. The differences in pointing caused by 
the atmospheric differential refraction are not impor- 
tant because the objects and the slit are relatively 
large. According to the slit size, the slit orientation, 
the airmass (1.09 < secz < 1.27), and the tables by 
Filippenko (1982), the difference between the sam- 
pled region by H/3 and that sampled by any other 
line is at most 8%; this difference will barely affect 
the line intensities because the object is considerably 
wider than the slit. 

We used a Th-Ar lamp for wavelength calibration 
in all spectral ranges and a tungsten bulb for internal 
flat-field images. The absolute flux calibration of all 
the spectra was achieved by taking cchcllograms of 
the standard stars HR 4963, HR 5501 and HR 7596. 
All of them are from the list of Hamuy et al. (1992), 
which includes bright stars with fluxes sampled at 
16 A steps. An average ciirvc for atmospheric ex- 
tinction was used (Schuster 1982). 

The spectra were reduced rising the IRAF^ 
echelle reduction package, following the standard 
procedure of bias subtraction, aperture extraction, 
flatflelding, wavelength calibration and flux calibra- 
tion. 

3. THE EMISSION-LINE SPECTRA 

The line fluxes and the dereddened intensities for 
NGC 5461 and NGC 5471, normalized to /(H/3) = 
100, are listed in Table 2. The line fluxes were mea- 
sured with the splot task of IRAF. The central 

''IRAF is distributed by NOAO, which is operated by 
AURA, under cooperative agreement with NSF. 
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OBSERVED AXD riEDDEXIXG-C^ORriEC^TED LINE RATIOS 









NGC 5461 






NGC 5471 




Ion 


Ao(A) 


Aobs(A) 


F(A)" 




Aote(A) 




lixr 


rn TTi 


O / Z U . V J o 


■3790 7n 


uz.uu 


oo.yu 


Q79Q ni 

/ zy .ui 


00.00 


A^ 01 
4:0. yi 


rn TTi 


Q790 QO 


0/0 1 


7fi 1 fi 

/ U. lU 


1 f)^ 40 
luu.^y 


Q7QI 7c: 


^1 91 

J- .Zl 


"^S ^4 

00. 0^4: 


TT T 

XI 1 


0000. uu 


0000. ZZ 


D.oy 


7 AS 


0000.00 


0.0 / 


f\ 9fi 

O.ZD 


ri\T„ TTTl 


ooOo. 1 


00 / l.OU 


10. io 


1 / .00 


00 / Z.ll 


fi^.OO 


AO ^A 
4y.04 


He T 
xlo 


QQQC f,^ ^ 
0000. UU 1 

oooy.uu J 




17.84 


23.48 


3892.20 


20.81 


23.22 


[IN C lilj 




oyTu.u^ 


fi9 
o.DZ 


A fi9 
4.DZ 


QQ7n OS 
oy ( u.yo 


1 /I HQ 
14. UO 


^^ A7 

iO.4 / 


H7 




oy 1 0.10 


1 ^ fi9 
10. OZ 


1 7 "^7 
1 ( .0 1 


^07^ i^S 

oy ( o.oo 


1 A n^i 

14. UO 


1 ^ AO 

10. 4y 


xie 1 




/inoQ t^n 
4uzy.ou 


l.Do 


Z.IZ 


AC\OQ S7 

4:Uzy.o 1 


9 O/l 
Z.U4 


9 OA 
Z.Z4 


[c; TTl 


/infis fin 


/in7i 7/1 


1 .zo . 


1 .00. 


/in71 07 

4u / 1 .y ( 


n 70- 
u. / y . 


n sfi- 

U.oO. 


N TTl 




/in70 "^n 
4u ( y.ou 


U.OO. 


U.OO. 


/1 070 79 


n fi/i • 

U.04. 


n 70- 

U. ( u. 


no 


/iini 7/1 

41U1 . 1 4 


Ai n'^ n^ 

'ilUu.Uo 


9"^ 

ZO. oy 


90 no 
zy .uy 


A1 n"^ AR 


9A A7 
Z4.4 / 


9fi fil 
ZO.Ol 






A'XAA 1 1 


A7 

■4:/ .00 


c^A f\7 


A'^AA CiS 

4044. uo 


Afi A"^ 
4D.40 


AO 1 9 
4y. iz 


rn TTTl 


40D0.UU 


/I'^fifi 7R 


1 r\A 


110 

i.iy 


40D ( .00 


S Sfi 
o.oO 


y .00 


Ho T 
lie i 


zi^sisi nn 


'ioyi.oo 


U.OO. 


n ^s- 

U.Oo. 








TTo T 
xlc 1 


/ l.^O 


AA7^ 9fi 


A A7 


A OS 


AA7^ 7(\ 
44 / 0. 1 D 


o.yo 


4 1^ 
4.10 


Hfi TT 


^OoO.uo 








Afion /ip; 

40yU.40 


n ^s 

U.Oo 


n "^0 
u.oy 


r A *. T"\/i 
[-rt.r 1 V J 










A71 ^ Sfi 
4 t 10. oO 


1 A^ 
1.40 


1 A^ 

1.40 


[-f\r 1 V J 


^ 1 ^u.zu 








A7A^ 1 n 

4 f 40. lU 


117 
1.1 / 


lis 

1. lo 


H 


/i^fii '^'^ 




1 nn nn 

lUU.UU 


1 nn nn 

lUU.UU 


/isfifi n^i 

4o00.U0 


1 nn nn 

lUU.UU 


1 nn nn 

lUU.UU 


Hr» 1 
nc 1 


/1 091 0'^ 


/109f^ 07 


1 Pin 

l.OU 


1 /IS 

1 .41:0 


/109fi Sf^ 

4yzo.oo 


n of^ 
u.yo 


n ot^ 
u.yo 


rn TTTl 


^yoo.yz 


/lQfi'4 90 


1 1 /I f^9 
1 14. OZ 


119 17 
1 IZ. 1 ( 


AOfi'^ sn 
4yoo.ou 


9'^9 Sfi 
ZOZ.oU 


9'^n OA 

zou.y4 


rn TTTl 


OUUO.oO 


p;ni 1 1 Pi 

OUll. iO 


000.01 


^9^ OS 

ozo. yo 


pini 1 7n 

OUll. / u 


(^SA 1 n 

Oo4.1U 


f>7f> 99 
D ( D.ZZ 


mi TTTl 
[v^l lllj 


uul ( .'^0 


OOZZ.UZ 


n ■ 

U.04:. 




^^00 09 
oozz . yz 


U.OO. 


n "^1 • 
u.oi. 


riv TTl 

[IN lij 


^7 


/ oy .00 


n AS 


n An 
u.^u 


( oy.4o 


U.OO. 


n ^9- 
u.oz. 


Ho T 


Oo 1 0.0 / 


OooU.OO 


1 fi9 
10. OZ 


1 ^ ^7 
10.0 / 


c^SSI 9S 
Oool.Zo 


IZ.OO 


1 1 ^9 

1 1 .OZ 


rn Ti 


R'^nn "^9 

OOUU. OZ 


DoUu. 11 


E^7 
0.0 ( 


/I 1 ^ 
4.10 


R'^OR 'iA 
D0UD.04 


/I /1 1 

4.41 


^ 09 

o.yz 


r^ TTTl 


fi'^i 9 nfi 

DOiZ.UO 


DOl / .00 


1 7R 
1. 1 


1 90 
i.zy 


R'^i 7 on 
001 1 .yu 


1 S9 
l.OZ 


1 fi9 
1 .OZ 


rn Ti 


0000. ol 




n 07 
u.y 1 


n 71 

U. / 1 


oooy.oy 


n on 
u.yu 


n sn 

U.oU 


r\T 1 ii 


fit^/is n'^ 

Du^o.Uo 


0000. oy 


91 n9 

Zl .uz 


1 /I on 
14. yu 








H^v 


DODz.oZ 


uOD / .00 


/ins '^1 


9SS ^7 
Zoo.O / 


0000.04 


■^91 S'^ 
OZl .00 


9sn 1 1 

ZoU. 1 1 


[Nil] 

L J 


6583.41 


6588.93 


59.33 


41.71 


6589.45 


8.62 


7.48 


He I 


6678.15 


6684.01 


4.80 


3.31 


6684.60 


3.73 


3.21 


[SII] 


6716.47 


6722.27 


29.42 


20.07 


6722.74 


12.70 


10.90 


[SII] 


6732.07 


6736.64 


22.90 


15.59 


6737.12 


9.44 


8.09 


He I 


7065.28 


7070.84 


3.22 


2.13 


7071.71 


2.85 


2.41 


[Ar III] 


7135.93 


7141.91 


15.52 


10.12 


7142.49 


7.91 


6.68 


[Oil] 


7319.65 


7326.60 


3.49 


2.22 


7327.25 


1.99 


1.66 


[Oil] 


7330.16 


7337.01 


3.03 


1.93 


7337.26 


1.96 


1.63 


[SlII] 


9068.90 


9076.42 


53.34 


28.47 


9077.07 


24.23 


18.84 


logF(H/3)'^ 




-12.41 
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0.18 ± 0.06 





''Normalized to H/3 = 100.00. ''Units of erg sec"^ cm"^. "Units of A. 
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wavelengths were measured with splot and referred 
to the hehocentric reference frame. The observa- 
tional uncertainties associated with the line flux in- 
tensities arc estimated to be 0.02 dex for those lines 
with F(A)/F(H/3) > 0.10, 0.04 dex for those lines 
with 0.02 < F{X)/F{HI3) < 0.10, and 0.08 dex for 
those lines with F(A)/F(H/3) < 0.02. Those lines 
marked with a colon are affected by an uncertainty 
of about 0.15 dex. 

Since the NIR and the blue spectra have been 
taken with different apertures, it is advisable to cor- 
rect for this effect (as well as for possible shifts in 
the telescope position from one night to another) 
to compute line ratios in a homogeneous manner. 
For this reason, we applied a grey shift to the 
line intensities of the NIR spectra lines, by match- 
ing the Ha intensities. The intensity of the blue 
spectrum has been considered for those lines falling 
in both wavelength ranges ([Ol] AA 6300, 6363, 
[Nil] AA 6548, 6584, [SIII] A 6312, and Ha). The 
shift applied was -F0.032 dex for NGC 5461 and 
-0.109 dex for NGC 5471. Note that for NGC 5461 
J(Ha)biue/-f(Ha)NiR > 1, even though the slit used 
in the blue range was smaller than the one used in 
the NIR: this fact probably implies that the telescope 
was not pointing exactly to the same position during 
the two nights, an effect amplified by the extremely- 
peaked brightness distribution of the region. Indeed, 
Figure 13 in Castafieda et al. (1992) and Figure 3 in 
Luridiana & Peimbert (2001) show that the Ha in- 
tensity falls to less than 10% of the peak value within 
3" from the center. 

The NIR lines were corrected for telluric absorp- 
tion; this correction was smaller than 3% for all the 
lines with the exception of AA7320, 9069 in NGC 
5461 where the correction amounted to 6% and 7%, 
and of A9069 in NGC 5471 where the correction 
amounted to 9%. The reddening coefficient, C(H/3), 
was determined comparing the /(H/3)//(Ha)biue ob- 
served ratio to the case B one computed by Storey 
& Hummer (1995) adopting = 8700 K and 
A/e = 150 cm-'^ for NGC 5461, and = 12,900 K, 
= 90 cm-3 for NGC 5471. The uncertainty 
in C(H/?) has been estimated assuming a 5% un- 
certainty in /(H/3) and /(Ha)biue- The results are 
shown in Table 2, together with the total uncorrected 
H/3 fluxes, measured in erg s~^ cm~^, the equivalent 
widths of H/3, and the reddening coeflicients. 

4. PHYSICAL CONDITIONS AND CHEMICAL 
ABUNDANCES 

The physical conditions in the two nebulae have 
been calculated with the ntplot and the temden 



IRAF (version 2.11.3) tasks. These tasks are based 
on 5-level model atoms for [Oil], [SIII], [CI III], and 
[ArlV], 6-level model atoms for [OlII] and [Nil], 
and an 8-level model atom for [SII]. The ntplot 
task allows a self-consistent computation of and 
Tg for a given ion, when the relevant line ratios are 
known; it was used to determine iVc(OlI), Te(OlI), 
and iVe(SlI), re(SlI) for both HII regions. When 
only one line ratio is known for a given ion, the 
task temden is used, which computes To for an as- 
sumed A^e value or vice versa. Te(NlI), re(OlII), 
and Te(SlII) were calculated with temden assum- 
ing Nc = A^c(OlI), and 7Vo(ArIV) was calculated 
assuming = re(OlII). It is worthwhile to re- 
mark that the high resolution and wide wavelength 
range of these observations considerably extends the 
number of available physical-condition diagnostics of 
the two H II regions, adding Ne{0 II), N^iAi IV) and 
To(S III) to the diagnostics studied in previous works. 

Table 3 lists the physical conditions of the two 
nebulae. The results obtained are in good agreement 
with those by Torres-Peimbert et al. (1989) within 
the estimated errors. The only notable exception 
is Te(SlI), for which Torres-Peimbert et al. (1989) 
found much higher values in both objects. However, 
such differences depend on the atomic parameters 
used, rather than on differences in the measured line 
ratios: indeed, using the intensity values of Torres- 
Peimbert et al. (1989) and the same atomic param- 
eters as us, one would obtain slightly lower tem- 
peratures than ours (the atomic parameters used in 
this paper are those from Verner, Verner, & Ferland 
1996, Keenan et al. 1993, and Ramsbottom, Bell, & 
Stafford 1996). 

The root mean square density, A^e(rms), to a very 
good approximation is given by the expression: 

_ f J(H/3) [l + N{Kc+)/N(ll+)]4nd^ \'^' 
Ae(rms)-| Vamn) hu^, j ' 

(1) 

where V is the observed volume, a(H/?,To) is the 
H/3 effective recombination coeflScient, and d is the 
distance to the region. We assume V = Anr^/S, 
where r is the radius of a circle equivalent in area 
to our slit, at a distance of 7.4 Mpc (Sandage & 
Tammann 1976, with the suggested correction by de 
Vaucouleurs 1978). The results, substantially equiv- 
alent to those of Torres-Peimbert et al. (1989), are 
also listed in Table 3. Prom these values, combined 
with /Vo(OlI) or iVo(SlI), we can calculate the fill- 
ing factor, which is 0.009 < e < 0.011 for NGC 5461 
and 0.013 < e < 0.022 forTvfG^C 5471. These values 
can be compared to the filling factor of the com- 
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TABLE 3 

ELECTRON TEMPERATURES AND DENSITIES OF NGC 5461 AND NGC 5471 



Quantity* 


NGC 5461 


NGC 5471 


Te(Oll) 

Tc(SlI) 

Te(NIl) 
Je(0 Illj 
Te(SlIl) 


10,400 ± 500 
8000 ± 1600 
8700 ± 700 
8500 ± 500 
9000 ± 800 


9500 ± 900 
9500 ± 1100 
8500 ± 550 
9300 ± 250 


14,200 ± 900 
10,000 ±3200 
10,800 ±2300 
13,000 ± 500 
12,100 ±800 


13,100 ± 2000 
12,700 ±2400 
10,800 ±1700 
13,400 ± 250 


A^c(OlI) 
A^c(SlI) 
A'"e(Ar IV) 
A'e(rms) 


150 ± 60 
130 ± 90 

14.2 


234 
14.8 


90 ±70 
70 ± ^[jo 
1350 ±^^0 
10.4 


186 
10.6 


References 


1 


2 


1 


2 



''Given in K and cm ^, respectively. References: — (1) This work; (2) Torres-Peimbert et al. 1989. 



positc model of NGC 5461 by Luridiana & Peimbert 
(2001), which averages 0.003. The lower filling factor 
derived by Luridiana & Peimbert (2001) was based 
on the data by Torres-Pcimbert et al. (1989), which 
obtained a higher A^c(S II) than that derived here. 

Small filling factors of ^ 0.01 arc characteristic of 
GEHRs (Shields 1986). Similar values of about 0.02 
have been found in the ionized gas of emission-line 
dwarf galaxies (Martin 1997). In contrast, Copetti 
et al. (2000) find a mean filling factor of 0.1 for a 
sample of 15 Galactic "normal" HII regions. Al- 
though the estimations of e in ionized nebulae are 
indirect and very rough, the difference between "nor- 
mal" and giant HII regions — if real — would imply 
a higher dumpiness in the giant ones. This is not 
unexpected considering the typical highly disturbed 
morphology and complex kinematics of GEHRs, that 
may be produced by the dynamical effects of stellar 
winds and supernovae in GEHRs. 

4.1. Ionic Abundances 
4.1.1. Heavy Elements 

The ionic abundances of heavy elements in the 
two regions have been calculated with the IRAF 
task abund, which is based on a three-zone model 
of the nebula, with the three zones corresponding to 
the low-, medium-, and high- ionization regions, each 
characterized by representative and values. 

Concerning the temperature, each of the 
three low-ionization diagnostics available — ^Te(NlI), 
To(SlI), and Tc(Oll) is affected by uncertainties: 
(a) Nitrogen lines are weak in these objects, es- 
pecially in NGC 5471. (b) Te(0 II) depends non- 
negligibly on the assumed reddening, and on the as- 



sumed density; furthermore, [O II] AA 7320, 7330 are 
weak and affected by telluric absorption lines. Fi- 
nally, (c) [S II] AA 4068, 4076 are also weak, and the 
estimates of the [S II] atomic parameters are con- 
stantly changing; additionally, To(SlI) has a depen- 
dence on the assumed N^. similar to that of rc(0 II). 
To minimize these uncertainties, we decided to adopt 
the average temperature of the three diagnostics, i.e., 
Te = 9000 K for NGC 5461 and = 11,700 K for 
NGC 5471. Finally, in both the medium- and high- 
ionization zones we adopted Te(0 III) as a represen- 
tative temperature in both regions, since re(S III) is 
more uncertain due to the observational errors in the 
sulfur lines and the frequent revisions of the [S III] 
atomic parametcirs. 

Concerning the electronic density, we adopted in 
the three zones of both HII regions A^e = A^e(0 II). 
We preferred iVo(OlI) over 7Ve(SlI), due to the fol- 
lowing facts: (a) the [O II] lines are more intense than 
the [S II] lines, (b) the [O II] electronic density diag- 
nostic is slightly more sensitive than the [S II] one 
in the low-density regime, and (c) the [S II] atomic 
parameters are still highly uncertain. 

With these assumptions, we have obtained the 
abimdanccs of N+, 0°, 0+, 0++, Ne++, S+, S++, 
C1++, Ar++, and At^+ (Table 4). 

4.1.2. Helium, 

There are 8 lines of the He I recombination 
spectrum observed in NGC 5461 and 7 lines in 
NGC 5471. To calculate the He+ abundance, it is 
necessary to correct the line intensities for the two 
well-known effects of self-absorption and coUisional 
enhancement. The correction for self- absorption 
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TABLE 4 

IONIC ABUNDANCES^ IN NGC 5461 AND 
NGC 5471 



Ion 


NGC 5461 


NGC 5471 


0° 


7 no _(_ n 2(1 


fi 63 + 20 


0+ 


8.08 ±0.14 


7.32 ±0.14 


0++ 


8.32 ±0.07 


8.02 ±0.07 


N+ 


7.04 ±0.12 


6.01 ±0.13 


Ne++ 


7.59 ±0.11 


7.32 ±0.10 


S+ 


6.02 ±0.19 


5.47 ±0.20 


S++ 


7.00 ±0.10 


6.45 ±0.10 


C1++ 


4.93 ±0.16 


4.20 ±0.16 


Ar++ 


6.15 ±0.12 


5.55 ±0.14 


Ar3+ 




5.07 ±0.10 



Given as 12 + log(X+VH+). 



may be rather uncertain for A 3889 and A 7065; 
furthermore the observational errors in the hne in- 
tensities of AA 4026, 4388, and 4922 are relatively 
large. Consequently we decided to determine the 
He+ abundances based on the AA 4471, 5876 and 
6678 lines. 

The He+ abundance can be computed for each 
He I line by means of the relation: 

+ = = ^ "eff(H/3) /(A)R 
^ ~ H+ 4861 aeff(A) /(H/J)' ^ ' 

where aeff(A) is the effective recombination coeffi- 
cient of the considered transition, and /(A)^ is the 
observed intensity of the considered line, corrected 
for collisional enhancement and optical depth effects. 

The collisional contribution to the helium lines 
was estimated from Kingdon & Ferland (1995) and 
Banjamin, Skillman, & Smits (1999) assuming = 
8700 K and = 150 cm'^ for NGC 5461, and 
Te = 12,900 K and = 90 cm'^ for NGC 5471. 
We decided to adopt for the temperature a weighted 
average of T(0 III) as representative of the 0++ re- 
gion and an average of T(OII), T(S II) and r(NII) 
as representative of the 0+ region; since the emis- 
sivity is dominated by the 0++ region the adopted 
temperatures are closer to the r(0 III) values. The 
self absorption effects in the triplet lines were es- 
timated from the computations by Robbins (1968), 
adopting r(3889) = 1.0 and 0.2 for NGC 5461 and 
NGC 5471, respectively. The t(3889) values were 
estimated from the AA 3889 and 7065 line intensities 
and from the photoionization model of NGC 5461 
by Luridiana & Peimbert (2001) computed with 
Cloudy (Ferland 1996). The aeff(A) coefficients 



TABLE 5 
HELIUM ABUNDANCE VALUES 



Line 


NGC 5461 


NGC 5471 


4471 


0.0983 


0.0846 


5876 


0.0945 


0.0867 


6678 


0.0836 


0.0879 


Mean value 


0.0937 


0.0864 



at Te = 8700 K (NGC 5461) and = 12,900 K 
(NGC 5471) have been obtained by power-law in- 
terpolation in the appropriate temperature inter- 
val between the data provided by Smits (1996) for 
iVe = 100 cm~^. Given the low densities of the two 
HII regions, the correction for collisional enhance- 
ment is generally small; similarly, the correction for 
optical depth effects is very small because r(3889) is 
relatively small. Finally, the effective recombination 
coefficient for H/3 is based on the data by Brockle- 
hurst (1971) for A^e = 100 cm~^, assuming a power- 
law interpolation for = 8700 K (NGC 5461) and 
Te = 12,900 K (NGC 5471). 

Table 5 presents the resulting y+ values, together 
with the mean value. The mean value was obtained 
by weighting the different lines as the square root of 
their intensities. 

There are three minor effects that were not 
considered in the helium abundance determination: 

(a) the correction due to the temperature structure, 

(b) the correction due to the collisional excitation of 
the Balmer lines and (c) the correction for the un- 
derlying absorption of the helium lines. To estimate 
the first two effects observations of higher accuracy 
are required (e.g., Peimbert, Peimbert, & Ruiz 2000; 
Peimbert, Peimbert, & Luridiana 2002; Luridiana, 
Peimbert, & Peimbert 2002), the third effect is con- 
sidered by Esteban et al. (2002) for these objects. 

In NGC 5471, Hell A 4686 is also observed, im- 
plying the presence of He++. The abundance of 
He++ relative to H+ can be computed by means of 
the equation: 

He++ _ 4686 Qon (H/3) /(4686) 
H+ ^ 4861 acff(4686) /(H/3) ' ^' 

where the Q;eff(A) coefficients at = 12900 K have 
been obtained by a power-law interpolation in tem- 
perature from the data provided by Smits (1996) 
for A'e = 100 cm^"^. The result of this calculation 
is He++/H+ = 3.9 x 10~^, a completely negligible 
quantity from the point of view of the total abun- 
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dance. It is interesting to note that in our spec- 
tra the nebular He II A 4686 hne is superposed on a 
much broader stellar line, about 3.5 times more in- 
tense than the nebular line; the sum of the two com- 
ponents is roughly a factor of 1.7 higher than the 
7(A4686) reported by Torres-Peimbert et al. (1989), 
while the nebular line alone is about one third of 
it. These differences are probably a combination of 
aperture effects, and the much lower resolution used 
by Torres-Peimbert ct al. (1989). 

The broad component of A 4686 is due to WR 
stars; in addition we also observe stellar features of 
NV 4602, 4620, CIV 4658, and CIV 5808, whereas 
N III 4640 and C III 4650 are not detected (see Ta- 
ble 6). All the detected lines carry uncertainties of 
the order of 20%. The NV 4602, 4620 feature im- 
plies the presence of WN stars, while the two C IV 
features implies the presence of WC or WO stars; the 
low 4658/5808 intensity rcatio, compared to the ob- 
served ratios compiled by Schaerer & Vacca (1998), 
might suggest that the WO stars dominate and the 
WC stars arc not numerous or even absent. However, 
it is not possible to interpret quantitatively the ob- 
served WR features in terms of a distribution among 
the various WR spectral types, due to the uncertain- 
ties both in our observations and in the calibration 
of WR spectra. Nevertheless, we can give a rough 
estimate of the total number of WR stars, using aver- 
age values for the "WR bump" luminosities for WN 
and WC stars (Smith 1991; Schaerer & Vacca 1998): 
We obtain in this way 14 "equivalent" WC4 stars 
and 26 "equivalent" WN7 stars; the WC4 stars ac- 
count for all the 5808 A bump flux, and almost one 
half of the N III/V + C III/IV + He II blend around 
4650 A, and the WN7 stars for the rest of the flux 
in the N III/V -|- C III/IV + He II blend. On the other 
hand, the observed H/? flux corresponds to a num- 
ber of ionizing photons (5(H°) ~ 5 x 10^^ s~^ (as- 
suming 0.8 for the covering factor), or 450 equivalent 
07V stars (Vacca 1994) ; for a Salpeter's IMF with 
Mup = 120 Mq, this corresponds to a total number 
of about 530 O stars in the region (Schaerer & Vacca 
1998) . From these estimates we derive a WR/ O ratio 
of 40/530 ~ 0.075. 

4.2. Chemical Abundances 

The chemical abundances for both H II regions 
were calculated under the assumption of no temper- 
ature fluctuations, and without considering the frac- 
tion of heavy elements trapped in dust. The result- 
ing values are listed in Table 7, which also contains 
the gaseous abundances derived by Torres-Peimbert 
et al. (1989) under the assumption that = 0, and 



TABLE 6 
STELLAR LINES OF NGC 5471 



A (A) 


Ion 




i{xr 


i(A)b 


4604 + 4620 


NV 


5.6 


9.0 


5.9 


4658 


CIV 


4.9 


7.8 


5.1 


4686 


He II 


3.9 


6.3 


4.1 


5801 + 5812 


CIV 


4.5 


6.5 


4.3 



'■Units: 10"^^ erg cm-^ s'^ ''Units: 10+^'^ erg s'^ 

the Orion and solar chemical abundances compiled 
from the literature. For both objects the Y values de- 
rived by us arc higher than those derived by Torres- 
Peimbert et al. (1989): we prefer our values due to 
the higher resolution of our spectra. For NGC 5461 
the N/0 and Ne/0 values derived by us are higher 
than those of Torres-Peimbert et al. (1989): the dif- 
ference is mostly due to differences in the icf's used 
by the two groups. For NGC 5471 the N/0 and 
Ne/0 values derived by us are similar to those by 
Torres-Peimbert et al. (1989), but in this case the 
ic/'s used were almost the same. The gaseous values 
for the heavy species in Orion are taken from Este- 
ban et al. (1998), with the exception of argon that 
was taken from Pcimbert (1993); the solar values arc 
taken from Holweger (2001) (N and Ne), Grevesse & 
Sauval (1998) (S, CI, and Ar), and from Christensen- 
Dalsgaard (1998) (He), while the oxygen abmidancc 
is an average of the values cited in AUende-Prieto, 
Lambert, & Asplund (2001) and Holweger (2001). 

4.2.1. NGC 5461 

The computation of the chemical abundances 

from the observed ionic abundances requires a 
knowledge of the ionization structure of the model. 
In the case of NGC 5461, we can rely on the ion- 
ization structure predicted by the photoionization 
model of Luridiana & Peimbert (2001 and private 
communication). Radial and volumetric integrations 
of the ionic fractions arc presented in Table 8. The 
radial average for a given ion X* is defined as 



X 



Rs 



N^dR 




N^dR, (4) 



while the volume average is: 



= J N{X')NedV I I N{X)N^dV. (5) 



Prom these average ionic ratios, the ionization 
correction factors can be easily computed, e.g., for 
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TABLE 7 

GASEOUS CHEMICAL ABUNDANCES IN NGC 5461, NGC 5471, ORION, AND THE SUN 





NGC 5461 


NGC 5471 


Orion 


Sun 


Element 




= 0) 


{-t 


= 0) 


(t^ = 0.024) 




Y 


0.285 ±0.011 


0.264 ± 0.007 


0.257 ±0.010 


0.241 ±0.007 


0.276 ±0.007 


0.271 ±0.010 


12 + logO/H 


8.52 ±0.10 


8.39 ± 0.08 


8.10 ±0.08 


8.05 ± 0.05 


8.64 ± 0.06 


8.71 ± 0.05 


logN/O 


-0.64 ±0.12 


-1.13 ±0.06 


-1.23 ±0.10 


-1.33 ±0.06 


-0.85 ±0.10 


-0.78 ±0.12 


log Nc/O 


-0.44 ±0.12 


-0.68 ±0.06 


-0.73 ±0.10 


-0.63 ±0.06 


-0.74 ±0.12 


-0.71 ±0.09 


logS/0 


-1.46 ±0.11 


-1.69 ±0.08 


-1.57 ±0.10 


-1.67 ±0.08 


-1.46 ±0.12 


-1.38 ±0.12 


logCl/0 


-3.54 ± 0.20 




-3.67 ±0.20 


-3.30 ±0.16 


-3.21 ±0.03 


-3.21 ±0.30 


logAr/0 


-2.35 ±0.12 


-2.21 ± 0.08 


-2.38 ±0.15 


-2.25 ±0.08 


-2.15 ±0.21 


-2.31 ±0.08 


References 


1 


2 


1 


2 


3, 4 


5, 6, 7, 8 



References:— (1) This work; (2) Torres-Peimbert et al. 1989; (3) Esteban et al. 1998; (4) Peimbert 1993; (5) Allende- 
Prieto et al. 2001; (6) Holweger 2001; (7) Christensen-Dalsgaard 1998; (8) Grevesse & Sauval 1998. 



the radial icf's: 



(X/XQ _ (H+/H);, 
(H/H+) (XVX)« ' 



(6) 



and analogously for the volume-weighted ionization 
correction factors. 

Helium. In NGC 5461 a non-nogiigiblc fraction of 
the total helium is in neutral form. According to the 
Luridiana & Peimbert (2001) model, such fraction is 
bracketed by the radial and the volume averages: 



0.02 < 



He" 
lie" 



< 0.27. 



(7) 



To estimate the He°/He amount we will take a 
weighted average of the radial and volume average 
ionization structure presented in Table 7. The rela- 
tive weights are 0.84 for the radial value and 0.16 for 
the volume average, and were obtained by adjusting 
the observed S+/S++ ratio to the ratio predicted by 
the model; this provides a good approximation since 
the ionization potential of He l is roughly equal to 
that of S II. From these relative; weights we obtain 
that He^/He = 0.929, and combining this value with 
the mean value of ?/+ (see Table 5) we find: 



or 



5 = y+/o.929 



4 X He/H 



0.1009, 



(8) 



0.285, (9) 



1 + 4 X He/H ± 2 X 16 X 0/H 

where the assumption has been made that oxygen 
represents one half of all the mass in heavy elements, 
and 0.08 dex were added to the gaseous value for 
0/H derived below to take into account the fraction 
of oxygen trapped in dust grains. 



Oxygen. The oxygen abundance was calculated 
by adding the observed ratios A^(0+)/7V(H+) and 
Ar(0++)/Ar(H+): 



O _ 0+ ± 0++ 
H ~ H+ 



(10) 



O*^ is not included in the sum since it is expected 
that the fraction of neutral oxygen is very similar to 
that of neutral hydrogen due to the charge-exchange 
reaction 0+± H" ^ 0"± H+ (see Osterbrock 1989, 
and Table 8 of this paper). 

For the case of 0+ and 0++, the ionic ratios 
amount to 0.354 and 0.644 (radial averages), or 0.708 
and 0.267 (volume averages), while the correspond- 
ing observed values are 0+/O*°*|obs = 0.352 and 
0++/O*°*|obs = 0.612. The fact that the observed 
values are in much better agreement with the radial 
than with the volume averages, reflects the fact that 
the slits used for the observations sample only a small 
part of the region. For this reason, we will consider 
in the following only the radius-weighted ionization 
correction factors, and we will drop for simplicity the 
superscript R, i.e., icf{X')^^°^ = ic/(X*)LPOi.«. 

The neutral oxygen is underpredicted by the 
model, no matter which kind of average is consid- 
ered (but particularly in the case of the radial aver- 
age). This fact is a common finding in photoioniza- 
tion modeling (see, e.g., Stasinska & Schaerer 1999; 
Luridiana & Peimbert 2001), but from the point of 
view of the computation of the ionization correction 
factors, it has no consequences since the 0° fraction 
is practically equal to the H° fraction. 

Nitrogen. Only two nitrogen ionic species, N+ 
and N++ , are expected in this H II region, and only 
N+ is observed. The model of NGC 5461 gives for 
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TABLE 8 

IONIZATION FRACTIONS FROM THE PHOTOIONIZATION MODEL OF NGC 5461"^ 

Ionization Stage 



Log (Radial Average) 



Log (Volume Average) 



Element 


I 


II 


HI 


IV 


I 


II 


HI 


IV 


Hydrogen 


-2.738 


-0.001 






-1.610 


-0.011 






Helium 


-1.661 


-0.009 






-0.567 


-0.137 






Nitrogen 


-3.045 


-0.836 


-0.069 


-3.605 


-1.785 


-0.270 


-0.351 


-4.288 


Oxygen 


-2.822 


-0.452 


-0.191 




-1.584 


-0.150 


-0.574 




Neon 


-2.842 


-0.169 


-0.494 




-1.815 


-0.058 


-0.960 




Sulfur 


-5.322 


-1.276 


-0.042 


-1.420 


-4.247 


-0.597 


-0.132 


-2.069 


Chlorine 


-4.716 


-1.093 


-0.041 


-2.037 


-3.644 


-0.477 


-0.177 


-2.689 


Argon 


-3.585 


-1.382 


-0.023 


-1.996 


-2.283 


-0.490 


-0.175 


-2.657 



*Luridiana & Peimbert 2001. 



the ionization correction factor of N+: 



ic/(N+) 



LPOl 



6.85, 



(11) 



where the radial average has been considered (see the 
previous section). It is interesting to compare this 
value with the one obtained for a "hot" region ac- 
cording to the method developed by Mathis & Rosa 
(1991): 

icf(N+)^^^^ = 3.31. (12) 

yielding for the total nitrogen abundance a value 
0.32 dex smaller than the one computed relying 
on the ionization structure of the model by Lurid- 
iana & Peimbert (2001). This difference is proba- 
bly related — at least partially to the consideration 
of the radius-averaged icf , which is less precise in 
the case of elements measured by means of lines ob- 
served with the longer slit. This hypothesis seems 
confirmed when one considers that the Mathis & 
Rosa (1991) result is bracketed by the radial and 
volume-averaged ic/'s. 

Neon. The Luridiana & Peimbert (2001) model of 
NGC 5461 gives: 

ic/(Ne++)LP°i = 3.12. (13) 



This value can be compared with the one pre- 
dicted by the method of Mathis & Rosa (1991): 



ic/(Ne++)^^9i = 1.64, 



(14) 



which gives again a total neon abundance about a 
factor of 2 smaller than the one based on the Lurid- 
iana & Peimbert (2001) model. In this case, the dif- 
ference cannot be explained in geometrical terms, as 



in the case of nitrogen and sulfur, since the volume- 
averaged ic/(Ne"''^) of the model by Luridiana & 
Peimbert, which amounts to 9.12. oven further away 
from the Mathis & Rosa (1991) value. This differ- 
ence should be studied further. 

Sulfur. In the case of sulfur, a correction must be 
applied to the ionic abundances of S"*" and S++ to 
account for the presence of a small fraction of S^~^ 
inside the H+ sphere. The Luridiana & Peimbert 
(2001) model gives: 



icf{S+ + S++)LPOi = 1.04, 



(15) 



yielding the total sulfur abundance listed in Table 7. 
This icf is almost equal to 1, implying that the S/H 
ratio is robust. In what follows, we explore the ad- 
ditional information for the ionization structure pro- 
vided by the S+/S ratio. 

The method by Mathis & Rosa (1991) gives: 



/(S+) 



MR91 



= 12.52, 



(16) 



while the photoionization model by Luridiana & Pe- 
imbert (2001) predicts 



k/(S+)'^™i = 18.88. 



(17) 



If the sulfur abundance were calculated from the ob- 
served S^ ratio only, applying the method by Mathis 
& Rosa (1991), the value logS/0 = -1.40 would 
be obtained. This value is 0.07 dex higher than 
the value derived before, obtained by direct obser- 
vations of both S+ and S^"*", combined with the 
icf{S++ S++) by Luridiana & Peimbert (2001). On 
the other hand, the value obtained by correcting the 
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observed S+/H+ by means of the jc/(S+) of Lurid- 
iana & Peimbert (2001) would be logS/0 = -1.22, 
a value 0.18 dex higher than the one in Table 7. An 
equivalent way of looking at the same discrepancy 
is to consider the S+ZS"*""*" value, which is 0.100 ac- 
cording to observations, and 0.053 according to the 
Luridiana & Peimbert (2001) model. This discrep- 
ancy can be qualitatively explained when one consid- 
ers that the S+ abundance value is based on obser- 
vations made with a larger slit than the one used in 
the case of, e.g., [Oil] A 3727 and [O III] A 5007 (see 
Luridiana & Peimbert 2001): the radial average is 
thus a better approximation in the oxygen than in 
the sulfur case. This view is supported by consid- 
ering that the observed S+/S++ ratio is bracketed 
by the radial average (S+/S"'""''!/? = 0.053) and the 



volume average (S+/S 



++\ 



0.253). However, as 



we have seen before, the icf's predicted by Cloudy 
are in general discrepant from those predicted by 
Mathis & Rosa (1991), so that also in this case fur- 
ther investigation is needed to assess the origin of 
the discrepancy. 

Chlorine. We adopted the icf of the model by 
Luridiana & Peimbert (2001) to determine the to- 
tal chlorine abundance, which is: 



-(-+\LP01 



1.10, 



(18) 



yielding a total chlorine abundance logCl/0 = 
—3.54. On the other hand, the value predicted by 
the method of Mathis & Rosa (1991) is: 



«c/(Cl++) 



MR91 



2.22, 



(19) 



which yields logCl/0 = —3.24, a value very close to 
the solar one (see Table 7). 

Argon. According to the Luridiana & Peimbert 

(2001) model, Ar++ accoimts for 95% of the to- 
tal argon, giving ic/(Ar++)^^°^ = 1.05 and im- 
plying the total abundance listed in Table 7, 
whereas the method by Mathis & Rosa (1991) gives 
icf{Ai++)^^^^ = 1.87, which would imply for the 
total argon abundance a value 0.25 dex higher than 
the adopted value. 

4.3. NGC 5411 

For NGC 5471 we do not have a detailed pho- 
toionization model available, thus we shall rely on a 
photoionization model of NGC 2363.^ NGC 2363 is 
a low-metallicity, high-excitation giant HII region, 

^More precisely, the model used is one of a series of simple 
burst models which constitute the extended-burst, composite 
model for NGC 2363. 



with a degree of ionization for the oxygen similar to 
that of NGC 5471. The model assumes Z = 0.20 Zq 
and a f = 3.0 Myr instantaneous burst, which ionizes 
a spherical region, made up of two concentric shells 
of different densities; see Luridiana et al. (1999) for 
further details. To improve the approximation, we 
did not directly use the ic/'s obtained from this 
model, but rather rescaled them to the actual ob- 
served ionic ratios. This procedure, which will be 
described more explicitly in the following sections, 
allows us to overcome possible differences related to 
the temperature and ionization structure of the two 
H II regions. The chemical abundances inferred for 
NGC 5471 are listed in Table 7, and will be discussed 
in detail in the following. 

Helium. The amount of neutral helium is prob- 
ably negligible inside NGC 5471, as implied by its 
high ionization degree (its ionization degree is simi- 
lar to that of M17 where it is found that the amount 
of neutral helium is negligible, see Peimbert, Torres- 
Peimbert, & Ruiz 1992; Estebanet al. 1999). There- 
fore, we find for the total helium abundance (see Ta- 
ble 5), 

He He+-hHe++ „ 

= 0.0868, (20) 



H 



or 



Y 



H+ 



4 X He/H 



l + 4x He/H 2 x 16 x 0/H 



0.257 (21) 



where the assumption has been made that oxygen 
represents one half of all the mass in heavy elements, 
and 0.08 dex were added to the gaseous value for 
0/H derived below to take into account the fraction 
of oxygen trapped in dust grains. 

Oxygen. The oxygen abundance was calculated, 
as in the case of NGC 5461, by simply adding the 

observed ratios 0+/H+ and 0++/H+, so no assump- 
tions on the ionization structure are necessary. 

Nitrogen. We determined the nitrogen abundance 
using the observed N+/0+ ratio and the 0+/0, 
N+/N ratios predicted by the model by Luridiana 
et al. (1999): 



o 



0+ 



obs 



O N+ / 



(22) 



LPL99 



Basically, this relation takes advantage of the simi- 
larity in the ionization structure of the two regions, 
rescaling to the observed 0+/0 for a better approx- 
imation. 

This equation yields for the total nitrogen abun- 
dance relative to oxygen the value log N/0 = 
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— 1.23 dcx. It is interesting to compare this value 
to the one obtained by means of the Mathis & Rosa 
(1991) method. The criterion of Mathis & Rosa 
places NGC 5471 into the category of "hot" H II re- 
gions, so that: 



icf{N+) 



MR91 



5.67, 



(23) 



which, combined with the observed value of N+ gives 
log N/0 = -1.34. On the other hand, with the 
usual assumption N+/0+ = N/0 one would obtain 
logN/0 = —1.31 dex. The difference between these 
two values and the one adopted by us is probably 
related to the non-negligible fraction of N"*"^ in the 
nebula, which is implicitly neglected by the empirical 
equation. 

Neon. In the case of Ne++, the scaling equation 
is: 



Ne 
O" 



Nc 



++ 



0^ 



/0++ Ne 



obs 



V O Ne++ 



(24) 



LPL99 



since the Ne"*""*" coexists with 0+"*". We obtain for the 
total neon abundance logNe/0 = —0.73 dex. This 
result is intermediate between the prediction of the 
Mathis & Rosa (1991) method, giving: 



k/(Ne++) 



++\MR91 



1.05, 



(25) 



and logNe/0 = —0.76, and the usual empirical 
relation Ne/0 = Ne++/0++, yielding logNe/0 = 
-0.70. 

Sulfur. The total sulfur abundance can be ob- 
tained from the following relations: 



o 



s++\ 

r 



0+ 



obs 



0^ 



O S++ 



and 



o 



V0++ 



/0++ s 



obs 



V o S++ 



LPL99 



LPL99 



(26) 



(27) 



since S++ coexists with 0+ and 0++. We adopted 
the value which is obtained by the weighted average 
of equations (26) and (27), where the weight is de- 
termined by the observed O ionization fraction, and 
obtained logS/0 = -1.43. 

The Mathis & Rosa (1991) method gives: 



/(S+) 



MR91 



= 17.68, 



(28) 



implying logS/0 = —1.38, in good agreement with 
the value derived with our model. 



Chlorine. We determined the chlorine abimdance 
with relations analogous to the ones used for sulfur: 



CI 
O" 



cr 



0+ 



X 
obs 



0+ _ci_\ 
"o" C1++; 



LPL99 



and 



CI 

o" 



\0++ 



X 
obs 



o 



++ 



CI 
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(29) 



(30) 
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and using the same weighting procedure used for sul- 
fur obtained logCl/0 = -3.71. The Mathis & Rosa 
(1991) method gives: 



*c/(Cl 

yielding logCl/0 = 
our adopted value. 



-I-I-\MR91 



12.35, 



(31) 



-2.81, a much higher value than 



Argon. A lower limit to the total argon abundance 
is given by: logAr/0 > log (Ar++-|- Ar^+/O)obs = 
—2.43. However, a correction should be made in or- 
der to account for the presence of a small fraction 
of Ar+. The equation based on the photoionization 
model gives: 
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^Ar+-FAr++ 
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0^ 
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O Ar+-hAr++ 
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(32) 

yielding logAr/0 = —2.38. This number can be 
compared to the prediction based on the Mathis & 
Rosa (1991) method: 



zc/(Ar 



(^) = 7-62, (33) 
\ -^r / MR91 



implying logAr/0 = —1.67, which is very large, 
probably meaning that the Ar++ fraction for this 
object is underestimated by the method of Mathis & 
Rosa (1991). 

5. DISCUSSION AND CONCLUSIONS 

We have described new high-resolution spectro- 
scopic data obtained for the two giant H II regions 
NGC 5461 and NGC 5471 in MlOl. We have 
determined five independent temperatures and four 

independent densities, as well as the abimdances for 
the following elements: H, He, N, O, Ne, S, CI, 
and Ar. 

The densities obtained from the forbidden lines 
are considerably higher than the root mean square 
densities, implying strong density fluctuations with 
typical values of 0.01 for the flUing factor. 
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The [ArlV] density for NGC 5471 seems to be 
larger than A'^e(OlI); since the [ArlV] Unes origi- 
nate close to the main ionizing sources, this find- 
ing probably indicates that the main ionizing stars 
are still located in a higher than average density re- 
gion. A similar result was obtained for NGC 2363 
by Perez, Gonzalcz-Delgado, & Vilchez (2001) and 
Esteban et al. (2002). 

A detailed comparison of the observed line in- 
tensities with those predicted by the photoionization 
model of Luridiana & Peimbert (2001) for NGC 5461 
was made. We found that the observed [O i] line 
intensities are higher than predicted by the model, 
indicating that the object is ionization bounded. A 
similar argument can be made based on the [S II] 
lines. 

The ic/(N+) for NGC 5461 derived from the pho- 
toionization model of Luridiana & Peimbert (2001) 
is about a factor of two higher than that derived by 
Mathis & Rosa (1991) and also about a factor of two 
higher than that derived from the empirical relation- 
ship given by N+ZO"*" = N/0. Further analysis is 
needed to understand the origin of this discrepancy. 

The zc/(Ne++) for NGC 5461 derived from the 
photoionization model of Luridiana & Peimbert 
(2001) is about a factor of two higher than that 
derived by Mathis & Rosa (1991) and also about 
a factor of two higher than that derived from the 
empirical relationship given by Ne^^/0+"'" = Ne/0 
(Peimbert & Costero 1969). In this case we think 
that the problem lies with the model by Luridiana & 
Peimbert (2001) for the following reasons: (a) prob- 
ably the best determination for the Ne/0 ratio of 
any HII region is that derived by Peimbert et al. 
(1992) for M17, amounting to -0.71 dex, (b) we 
expect the Ne/0 abundance to be practically the 
same for NGC 5461, NGC 5471, and M17 because, 
to a very good first approximation, these elements 
arc formed by massive stars and arc ejected into the 
interstellar medium by supernovac of type II, and 
(c) the jc/(Ne++) for NGC 5471 and M17 are con- 
siderably smaller than for NGC 5461 indicating that 
the abundances for the first two objects arc more 
reliable than for NGC 5461. The large Ne/0 value 
derived for NGC 5461 is probably an upper limit 
due to the presence of the charge-exchange reaction 
0+++HO 0++H+ that allows some 0+ to coex- 
ist with Nc++; this effect will reduce the zc/(Ne^^) 
predicted by the model (see Peimbert, Luridiana, & 
Torres-Peimbert 1995, and references therein). 

For NGC 5471, the ionization correction factors 
were computed based on a photoionization model 
for a similar HII region. Comparison of the total 



chemical abundances obtained to those implied by 
the Mathis & Rosa (1991) method highlighted sensi- 
tive discrepancies, particularly in the case of chlorine 
and argon. The origin of these discrepancies requires 
further study. 
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